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We study the interplay between strong correlations and in- 
commensurability on fermions using mean field as well as ex- 
act many-body Lanczos diagonalization techniques. In a two- 
dimensional parameter space, mean field phase diagram of 
infinite system shows a critical phase with multifractal char- 
acteristics sandwiched between a Bloch-type extended and a 
localized phase. Exact numerical computation on finite size 
systems of the Kohn charge stiffness D c , characterizing trans- 
port properties, and the charge exponent K p , characterizing 
the superconducting pairing fluctuations, shows the existence 
of a phase with superconducting correlations. This phase may 
be characterized by a power law scaling of the charge stiffness 
constant in contrast to the Anderson localized phase where 
D c scales exponentially with the size of the system. We argue 
that this intermediate phase may be a dressed up analog of 
the critical phase of the noninteracting fermions in an incom- 
mensurate potential. 



I. INTRODUCTION 

The study of metal-insulator transition in strongly cor- 
related one-dimensional (ID) disordered systems have 
been the subject of many recent studies. HQ These stud- 
ies have been motivated by various different viewpoints. 
A great deal of theoretical research in this area has been 
to account for the magnitude of persistent currents in 
mesoscopic rings. |3||| Another motivation for these stud- 
ies has been the argument that ID models in a staggered 
field may describe the physics of two-dimensional sys- 
tems; the subject of great interest in high T c supercon- 
ductivity. || 

In this paper, we describe our study of the phase dia- 
gram of many-body spinless fermion chain in presence of 
a sinusoidal potential with periodicity which introduces 
a new length scale in the problem. We approach the 
problem of understanding the phase diagram of this ape- 
riodic system from a perspective very different from the 
previous studies. It has been known that noninteracting 



systems in incommensurate potentials exhibit metallic, 
insulating as well as an exotic fractal phase which is in- 
between metallic and insulating and is usually referred 
as critical. These systems which are in-between pe- 
riodic and random have attracted a great deal of atten- 
tion due to the possibility of exhibiting a novel type of 
transport as the multifractal wavefunctions are charac- 
terized by almost a power law decay. Technical complex- 
ities associated with studying strongly correlated systems 
make the study of the effect of strong correlations on in- 
commensurate systems a rather difficult task. Hence, 
the natural questions like how the metal-insulator tran- 
sition and the fractal characteristics of non-interacting 
fermions are affected when many-body fermion-fermion 
correlations are taken into account has remained an open 
theoretical challenge. 

We recently published our preliminary results Q 
on this problem using Lanczos diagonalization method 
where we focussed on the transport studies of ID spinlcss 
model in the presence of both repulsive as well as the at- 
tractive interaction. One of the intriguing results of our 
study for attractive interaction was the possibility of a 
new phase where the charge stiffness may scale as a power 
law. The purpose of this paper is three fold. Firstly, we 
present our simulation results for systems of bigger sizes 
than those reported in our earlier paper. This strength- 
ens the possibility of a new type of phase in strongly cor- 
related systems. Secondly, we compute the charge expo- 
nent K p which characterizes the long range correlations 
in the system. Our results show that the regime with 
superconducting fluctuations overlaps the regime where 
we see a new type of transport characteristics. Thirdly, 
we describe self-consistent mean field (MF) theory results 
where the mean field is assumed to be due to copper pairs. 
The self-consistent mean field theory shows the existence 
of a critical phase, sandwiched between an extended and 
localized phase. Interestingly, this critical phase exists in 
a finite parameter space of measure directly proportional 
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to the strength of the cooper pair interactions. We cor- 
relate the MF and the exact diagonalization results and 
argue that the new type of behavior seen in numerical 
simulation may be a dressed up analog of the critical 
phase of noninteracting fermions. 

In section II, we describe our model and two interesting 
limits where the model has been extensively studied. In 
section III, we discuss mean-field results for the spinless 
fcrmion case. In addition to the various remarks in the 
introduction, this section will further bring to focus our 
motivations for the proposed study. In section IV, we de- 
scribe our results for the charge exponent K p and charge 
stiffness D c . In section V we summarize our results. 



II. MODEL SYSTEM AND ITS TWO 
INTERESTING LIMITS 

Consider an interacting spinless fermion model on a 
ID ring in a quasiperiodic potential, 
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H = - ^(cjcj+i + cj +1 c 4 ) + V n, L n t+1 + ^ h, L n t . 

i—1 i—1 i—1 

(1) 

The Ci and are respectively the fermion and number 
operators at site i. The site dependent potential is chosen 
to be of the form, hi = A cos(27reri). Here, A represents 
the strength of the potential and er is an irrational num- 
ber which is chosen for convenience to be the Golden 
Mean ( V ^T 1 )■ This parameter describes the two com- 
peting length scales in the system: namely, the fermion 
lattice space and the periodicity of the sinusoidal poten- 
tial. 

The above model has two interesting limits that have 
been studied extensively. For A = 0, the interacting spin- 
less fermion problem could be mapped to the Heisenberg- 
Ising XXZ spin proble m fis| for which a closed Bethe's 
ansatz solution exists. |l9| For all densities, the model 
exhibits a Luttinger liquid phase. The Luttinger liquid 
phase differs from the conventional Fermi liquid which ac- 
counts for the conduction behavior in conventional metal- 
lic phase. In Fermi liquid phase the low energy excited 
states of the interacting electron gas can be described in 
term of quasiparticle excitations which are analog to the 
single particle excitations of a free Fermi system. In con- 
trast, the Luttinger phase is characterized by many-body 
collective excitations. For attractive interaction (V < 0), 
the model exhibit SC correlations that can be described 
in terms of the charge exponent K p , 



< OiO l+r >oc r 



where Oi is a (singlet) pairing operator, 



Oi = 



(cjC i+ i + h.c.) 



(2) 



(3) 



Therefore, K p > 1 signals the existence of SC fluctua- 
tions. 

Other interesting limit of the model is the V — case 
where the model describes the famous Harper equation 
p0| which has also been studied both numerically as well 
as more recently by Bethe's ansatz. The Harper equa- 
tion is a paradigm in the study of low-dimensional in- 
commensurate systems as it exhibits a metal-insulator 
transition in one dimension. At the onset of transi- 
tion A c = 2, the quantum states are neither extended nor 
localized but instead exhibit fractal character and have 
been termed as critical. The spectrum contains an infi- 
nite number of gaps and is believed to be a Cantor set 
of zero measure. These interesting aspects of the wave 
function and the spectra have been shown to be reflected 
in the transport properties such as Launder resistance. 
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The possibility of a phase which is neither Anderson 
insulating nor metallic in strongly correlated systems is 
a fascinating theoretical problem. We explore this first 
using mean-field theory and then using exact many body 
Lanczos diagonalization on finite size systems. In cor- 
related systems, we will focus directly on the transport 
properties characterized by the stiffness constant. This 
is in contrast to the mean-field case where the single par- 
ticle wavefunctions can be used as the diagnostic tool for 
characterizing the nature of the phase in the system. 



III. MEAN FIELD THEORY 

For attractive interaction, we assume the existence of 
electron pairing resulting in a mean field A 
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< CiC l+ l > 
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In terms of the mean field A, the Hamiltonian becomes 
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(5) 



The parameter A is a function of both t and V and has 
to be determined self-consistently. 

The above mean-field Hamiltonian is the fermion rep- 
resentation of spin-i XY chain in a transverse field. |23| ] 
The fermion operators c, are related to spin operators Si 
by Jordan Wigner transformation, 



Sn = c l exp(i7r ^2 c lnC m ) 
S~ = exp(-Z7T ^ 4n c m)c n 



(6) 
(7) 



2 



S: 



-A, 



(8) 



and for 

Z7T 



The parameters t and VA are related to the exchange 
interactions in the spin model: t is proportional to the 
sum of the exchange interaction in the 2-dimensional spin 
space while the fermion non-conserving term VA is pro- 
portional to spin space anisotropy. The V dependent 
term in the Hamiltonian breaks the 0(2) symmetry of the 
spin chain and the spin model is in the universality class 
of the Ising model with anisotropy proportional to VA. 
The mean-field A can be determined self-consistently by 
the methods described by Lieb et al (for example as the 
results shown in fig. 1). |23j In the absence of magnetic 
field (A = 0), the A can be determined analytically in 
the following two limits: For y = 0, A 
t = AV, A = |. 

The anisotropic spin model in the presence of modulat- 
ing magnetic field was studied recently as the perturbed 
Harper equation, jnj Numerical diagonalization as well 
as exact decimation scheme showed that the model ex- 
hibits extended (E), localized (L) as well as a critical 
(C) phase. A novel aspect of the model is the fact that 
unlike Harper, the C phase exists in a finite parameter 
interval. The term VA results in fattening the critical 
point of Harper to a critical phase, sandwiched between 
the Bloch phase and the localized phase. The boundaries 
between the Bloch states and critical states as well as the 
boundaries between the critical and localized states were 
found numerically to be determined by simple relations: 

• A < 2(t + VA) , we have a E phase. 

• If A > 2[t — VA) , we have a L phase. 

• If 2(t + VA) < A < 2(t + VA) , we have a C phase. 

Therefore, by determining the parameter A self- 
consistently, the phase boundaries of the model can be 
inferred. It should be noted that the mean-field theory 
can also be done in the absence of cooper pairs by writ- 
ing mean field term as < cjcj+i >. In this case, the MF 
Hamiltonian belongs to the universality class of Harper 
equation which exhibit E-L transition. The critical phase 
in this case reduces to a critical line. Such a MF theory 
gives a self-consistent solution for both repulsive as well 
as attractive interaction V. 

Figure 1 shows the mean field A obtained self- 
consistently. It should be noted that the strength of the 
aperiodic field A determines a critical value of V below 
which the mean field A is zero. For repulsive interac- 
tion, the self-consistent value of A was always found to 
be zero. 

Figure 2 shows the MF phase diagram in the two- 
dimensional parameter space of V — A. For attractive 
interaction, the model exhibits E, L as well as fat C phase 
beyond a critical value of A. Exact renormalization study 



of the C phase has shown that the wave function ex- 
hibits self-similarity. The exponents characterizing the 
self-similar behavior at the onset of E-C and C-L transi- 
tion are different from those of the fat critical phase sand- 
wiched between these two phases. |[(| Another interest- 
ing aspect of this diagram is the fact that C-L transition 
line coincides with the magnetic transition to long range 
order which is in the universality class of Ising model in 
transverse field. 

The MF model exhibits a gap in the energy spectrum 
(except at the conformal point which corresponds to the 
transition to magnetic long range order in the XY spin 
model). However, we argue that the transition under 
consideration in disordered systems is the Anderson tran- 
sition and therefore the transport properties are deter- 
mined by the nature of single particle wave function. 

The Lanczos's numerical results shown in the next sec- 
tion confirm the existence of superconducting pairing for 
attractive interaction. We would like to argue that this 
justifies our basic assumption about the existence of a 
mean field is derived from the cooper pairs. Reentrant 
nature of the MF diagram may be just the artifact of 
MF assumption. Therefore, the MF phase diagram may 
describe qualitatively the physics of strongly correlated 
systems in the regime where the superconducting fluctu- 
ations exist. In spite of the limitations of MFT in ID, 
the interesting phase diagram provides a strong incentive 
for further exploration of phase diagram with a view to 
check the possibility of an exotic phase (with localization 
properties intermediate between metallic and insulating) 
in strongly correlated systems. 



IV. LANCZOS DIAGONALIZATION RESULTS 

Motivated by the possibility of a new type of phase, we 
next use Lanczos exact diagonalization method to obtain 
the stiffness constant D c which determines the transport 
properties of the system. In addition, we also compute 
the charge exponent to check the validity of MF theory 
based on the assumption that the MF derived from the 
fermion pairing. As explained below, both these quanti- 
ties can be determined from the ground state properties 
of the system. 

The first step in this computation is to determine the 
ground state energy Eq of a one-dimensional ring with 
N lattice sites subjected to a transverse magnetic flux 
<F The the Kohn stiffness constant D c is then given by, 
§§ 



N rf 2 g ($) 



(9) 



The charge exponent K p can also be determined in terms 
of ground state properties by first computing the com- 
pressibility k. As explained by P. Prelovsek et all || the 
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compressibility for finite systems can be calculated from 
values of ground state energies: 

1 N 

— = — [E (N e + 2) - 2E (N e ) + E (N e - 2)] (10) 
p z n 4 

Here N e are the number of fermions in the system with 
fermion density p — N e /N. Now the following two equa- 
tions relate k to the the charge exponent via the the 
renormalized Fermi velocity u p , 

K p = 7T— — — (11) 



Our simulations are done for systems of various sizes 
N and electronic densities p = in two-dimensional 
parameter space V and A. To simulate golden mean 
quasiperiodicity into the model, we used Lanczos meth- 
ods for systems of various Fibonacci sizes. Furthermore, 
we worked with densities which are the rational approxi- 
mants to the the square of the golden mean a 2 . This pro- 
cedure provides several possible sizes (5, 8, 13 and 21) for 
which the Lanczos diagonalization can be done keeping 
the density of the fermions almost a constant. It should 
be noted that unlike the previous studies involving ran- 
dom disorder, we cannot work with arbitrary sizes. This 
limits not only the number of sizes that we can study, 
but also forces us to work with densities different from 
half-filling . Therefore, our studies are for systems away 
from half-filling where the umklapp processes become ir- 
relevant and the system in absence of disorder is metal- 
lic. The next possible Fibonacci size 34 is rather hard to 
simulate using present day technology, and therefore, our 
results are in fact for the maximum possible exact diag- 
onalization size (for simulating incommensurate effects) 
that can be done with current regular computers. To ob- 
tain additional data points, we also show some results for 
densities which are rational approximants of ^"g -1 such 
us p = y|. However, the fermion density for this case is 
different from those of Fibonacci size systems. 

Figures 3 and 4 respectively show our results for the 
charge exponent and the charge stiffness. Figures (a) and 
(b) correspond to two different densities and are included 
here to show the consistency of our results independent 
of the density. As shown in figure 3, comparison with 
the A = case show that the presence of a new compet- 
ing length enhances the SC pairing fluctuations in the 
strongly correlated fermions. This effect becomes spe- 
cially strong around V — — 2 where the charge stiffness 
attains a maximum value as shown in figures 4 and 5. 
The existence of a characteristic peak in D c was reported 
in our earlier paper for N < 13. Simulation results pre- 
sented here confirms our earlier results for bigger size 
systems. What is new here is the fact that by computing 
the charge exponent K p , we are now able to correlate 



the regime characterized by a clear peak with the regime 
where the model exhibits SC fluctuations. This also pro- 
vides a plausible argument to justify the use of an effec- 
tive mean field approach (like the one use in section III) 
for the Copper like quasiparticles. Since the MF phase di- 
agram exhibits a critical behavior with fractal self-similar 
wave function, we speculate that the phase diagram in 
the vicinity of the characteristic peak with strong SC 
fluctuations is a dressed up version of the critical phase 
in systems with strong fermion correlations. We would 
like to stress that we view our MF results as providing 
only qualitative predictions as the phase boundary (C-L) 
predicted by MF does not agree with the one obtain by 
the exact diagonalization method. 

Unlike the MF case where the intermediate nature of 
the phase was due to the fractal nature of the single parti- 
cle wave function, the many body phase diagram is char- 
acterized by its transport characteristics. Figures 4 and 5 
indicate the possible existence of a region where D c may 
take intermediate values: between those of a metallic and 
those corresponding to the Anderson localized insulating 
phase. Our simulations show that the height of the peak 
in D c decreases rather slowly with the size of the sys- 
tem. This is in contrast to the Anderson localized phase 
where the D c decays exponentially with the size of the 
system. We conjecture that in the regime near the peak, 
the charge stiffness decays as a power law. This conjec- 
ture was verified only at a at the special point V = 0. 
For arbitrary value of V, it is rather difficult to verify 
this conjecture for large N. 

V. CONCLUSIONS AND DISCUSSION 

Two central results of this paper are: the existence of 
SC fluctuations in the strongly correlated fermion model 
with competing length scales where the incommensura- 
bility enhances the SC fluctuation and a possibility of an 
intermediate phase (in-between metallic and Anderson 
localized) which may be related to the critical phase of 
nonintcracting model which exhibits fractal characteris- 
tics. This is the first paper where the effects of two in- 
commensurate lengths are studied in strongly correlated 
system. We hope that our results will stimulate further 
studies of this type of behavior in other models such as 
t — J and Hubbard models. 

Our numerical calculation of the charge exponent pa- 
rameter describing SC pairing coherence provides a jus- 
tification of our mean field ideas where fermion-fermion 
correlations were explicitly assumed. We want to empha- 
size that in aperiodic case, one cannot disregard the mean 
field results in ID as one does in the pure models. This 
is because, in the pure model, the metal-insulator tran- 
sition is the Mott transition where a conducting phase 
becomes insulating due to the opening of a gap. Previ- 
ous studies have shown that mean field theory contra- 
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diets the exact results as far as the existence of a gap 
is concerned. In the aperiodic cases, the metal-insulator 
transition is the Anderson transition. In the noninteract- 
ing model, this transition is characterized by the local- 
ization of the single-particle wave function and a point 
spectrum. Therefore, we would like to argue that even 
though mean field results are incorrect regarding the ex- 
istence or nonexistence of a gap, it may still be of some 
validity in describing the Anderson transition. 

Recently, there have been density matrix renormaliza- 
tion group (DMRG) studies of spinless fermion models 
that show the existence of a delocalized phase for Ander- 
son disordered potentials. (25) The range of the parame- 
ter V for which this phase seems to exist coincides with 
the values of V for which we postulated the possibility 
of an intermediate phase. We think that these findings 
further support our argument that an intermediate phase 
may exist for the case of quasiperiodic potentials. 

Our previous results, hinting a new mechanism in- 
volving some sort of competition other than the known 
screening of the disorder due to SC fluctuations, are 
strengthened by our studies with bigger size systems. 
Furthermore, the study, involving mean field and exact 
diagonalization, provides a more convincing argument 
that this new phase may be related to the critical phase 
of non-interacting systems. 
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FIG. 1. Self-consistent mean-field A vs V for several val- 
ues of A. Namely, solid line A = 0.5, dotted line A = 1, 
short-dashed line A = 1.5 and long-dashed line A = 2.1. 

FIG. 2. Mean-field phase diagram in V — A plane showing 
the E, the C and the L phases. Along the dashed line describ- 
ing the C-L transition, the model is conformally invariant. 
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FIG. 3. The Charge exponent K p vs V for A = 0.0 (dashed 
curve) and A = 1.0 (un-dashed curve). Part (a) for p = y| 
and part (b) for p = || showing that incommensurability 
enhances the SC pairing. 

FIG. 4. Charge stiffness D c versus V for A = 1.0. Part 
(a) shows p = part (b) p = || . Shaded parts indicate the 
regime where SC fluctuations exist. 

FIG. 5. Charge stiffness versus V for p — i| . and A = 1.0. 
The points obtained for this case are indicated by crosses. The 
interpolated dashed curved is mean to be used as a guide to 
the eye. 
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